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Abstract. 

An exciting unsolved problem in the study of high energy processes of early type stars concerns the physical 
mechanism for producing X-rays near the Be star 7 Cassiopeiae. By now we know that this source and several 
"7 Cas analogs" exhibit an unusual hard thermal X-ray spectrum, compared both to normal massive stars and 
the non-thermal emission of known Be/X-ray binaries. Also, its light curve is variable on almost all conceivable 
timescales. In this study we reanalyze a high dispersion spectrum obtained by Chandra in 2001 and combine it with 
the analysis of a new (2004) spectrum and light curve obtained by XMM-TVeuiton. We find that both spectra can be 
fit well with 3-4 optically thin, thermal components consisting of a hot component having a temperature kTg ~ 
12-14 keV, perhaps one with a value of ~ 2.4 keV, and two with well defined values near 0.6 keV and 0.11 keV. We 
argue that these components arise in discrete (almost monothermal) plasmas. Moreover, they cannot be produced 
within an integral gas structure or by the cooling of a dominant hot process. Consistent with earlier findings, we 
also find that the Fe abundance arising from K-shell ions is significantly subsolar and less than the Fe abundance 
from L-shell ions. We also find novel properties not present in the earlier Chandra spectrum, including a dramatic 
decrease in the local photoelectric absorption of soft X-rays, a decrease in the strength of the Fe and possibly of 
the Si K fluorescence features, underpredicted lines in two ions each of Ne and N (suggesting abundances that 
are ~ 1.5-3x and ~ 4x solar, respectively), and broadening of the strong NeXLya and OVIIILya lines. In 
addition, we note certain traits in the 7 Cas spectrum that are different from those of the fairly well studied 
analog HD 110432 - in this sense the stars have different "personalities." In particular, for 7 Cas the hot X-ray 
component remains nearly constant in temperature, and the photoelectric absorption of the X-ray plasmas can 
change dramatically. As found by previous investigators of 7 Cas, changes in flux, whether occurring slowly or in 
rapidly evolving flares, are only seldomly accompanied by variations in hardness. Moreover, the light curve can 
show a "periodicity" that is due to the presence of flux minima that recur semiregularly over a few hours, and 
which can appear again at different epochs. 
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' 1. Introduction 

• The observed properties of 7 Cas (B0.5 Ve; my = 2.25) in 
[ ultraviolet to infrared wavelengths have led to major dis- 
I coveries related to the Be-phenomenon since its discov- 
' ery as the fi rst of what became known as "Be stars" 
( Secchil 118671 ). However, its X-ray emission is not typ- 
ical of this kind of object (jWhite et al.lll982l ). The X- 
ray emission of 7 Cas is dominated by a hot plasma 
component {kT ^ 10-12 keV), with mildly high lumi- 



nosity ('^10'^^ "^^ergs ^) and variable flux (jSmith et al 



2004 and references therein). Be stars at large are softer 
X-ray emitters {kT ~ 0.5 keV) with l ower luminosity 



10 ^^ ergs ^) and little variability (e.g.. lBergh6fer et al 



19971 ). In contrast, most of the well investigated Be/X-ray 
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binary systems are Be-|-neutron star systems and show 
a high luminosity (^ 10'^'^ ergs~^) and a distinctly non- 
thermal high energy distribution. Notably, no large X-ray 
outburst has been observed in 7 Cas. This is at variance 
with the behavior witnessed in many classical Be/X-ray 
systems. 7 Cas is now known to be in a 204-day binary, 
with an ecc entricity variously dete r mined to be near e = 
and 0.26 ( Harmanec et al. 2000l : iMiroshnichenko et al. 
I2OO2I) . Little is known about the secondary, except that it 
is likely to have a mass in the range 0.5-2 Mq. Recently, X- 
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ray properties similar to those of 7 Cas have been observed 
in a s mall but grow i ng number of Be stars: the 7 Cas-like 



stars (iMotch et all l2007t iLooes de Oliveira et all l200d: 



Lopes de Oliveirall2007tr 

A conclusive explanation for the X-ray emission of 
7 Cas and its analogs is lacking. The first suggestion to 
explain the hard X-ray emission of 7 Cas was that it is 
powered by accretion of matter from the Be wind or disk 
onto an accreting neutron star (jWhite et al.lll982l ). More 
recently, some have advocated the presence of an accret- 
ing white dwarf companion. This suggestion comes from 
several X- ray characteristics of cataclysmic variab le sys- 
tems (e.g., Murakami et al. 1986t Kubo et al. 19981) . such 
as their thermal nature and only moderate X-ray lumi- 
nosity. However, the analogy is incomplete upon further 
scrutiny, and X-rays would have to come from accretion 
onto a white dwarf representative of a new c lass of Be/X- 



ray binaries (jLopes de Oliveira et al.l |2007|) and have 



high efficiency of mass accretion energy to X-ray flux. 

A second proposed scenario suggests that the X-ray 
emission of 7 Cas is a consequence of magnetic interac- 
tion between its stellar surface and a Be (Keplerian de- 
cretion) circumstellar disk that entrains a magnetic field. 
This idea is based in part on the correlations of the UV 
and optical variabilities with the X-ray light curve and 
also indirect evidence of magnetic field observed in this 
star. This evidence takes the form of migrating subfea- 
tures running blue-to-red through optical and UV line pro - 
files (|Yang et al.lll988t ISmith. Robinson, fc Hatzeslll998l) 
and also the discovery of a gray, robust, 1.21581-day 
feature in the star's light curve. This periodicity must 
be very near or identical to th e star's rotation period 
( Smith. Henry, fc Vishniac 2006f) . The inference of a disk 
connection is based on the cyclical aspect and reddish 
tinge of the optical cycles that correlate with the X-ray 
ones. It is also supported by patterns of variability in the 
C IV and Si IV lines that are occas ionally seen in the X-ray 
light curve ( Cranmer et al. 2000l ). One piece of evidence 



for this interaction comes from the observation of highly 
redshifted absorption UV lines, suggestive of material be- 
ing ejected from the circumstellar environment toward the 
star with enough energies to produce ~ lOkeV X -rays 
when they impact the star ( Smith fc Robinson 1999t ). 

Both the magnetic disk and accretion interpretations 
have important astrophysical implications. The former 
suggests that disk dynamos operate and also that these 
stars are possible proto-magnetars. The accretion inter- 
pretation would imply the presence of a neutron star in 
an unusual accretion regime, or a white dwarf with novel 
properties. The presence of a white dwarf is still specu- 
lative, but such objects are predicted to be common as 
secondaries in models of evolution of binaries with B pri- 
maries. The resolution of the mystery of the production 
of X-rays in "7 Cas stars" could lead to a breakthrough in 
either of these fields. 

We have focused our efforts to try to understand the 
origin of the X-rays of 7 Cas and analog Be stars. In this 
paper we confine ourselves to the X-ray properties of 7 Cas 



from high and medium resolution spectroscopy and tim- 
ing studies using data obtained by XMM-Newton satellite. 
Chandra HETG data are also reinvestigated. 

2. Previous X-ray observations of 7 Cas 

Several Rossi X-ray Timing Explorer (RXTE) 
observations of the light curve o f 7 Ca s de- 
Smith. Robinson, fc Cor"be^ (|l998l) and 
(hereafter "SRC98" and 



tailed bv 

Robinson fc SmithI (|200oF 



"RSOO," respectively) have disclosed that the light curve 
undergoes variations on rapid (flaring), intermediate 
(several hours) and long (2-3 month cycles) timescales. 
We summarize the flaring results first as follows: 

— Flares (shots) are ubiquitous, except during brief pe- 
riods of very low X-ray flux. The reoccurrence of these 
lulls is often cyclical. 

— Individual flare profiles are narrow and symmetrical in 
shape. 

— Collectively flares show a log-normal distribution in 
energy. 

— The flares show an approximately 1/f distribution 
down to the photon limit of the instrument (for RXTE 
about 4 seconds). Occasional groups of flares (aggre- 
gates) last as long as a few minutes, and indicate the 
end of this red noise dependence. 

— Spectrophotometry of shots indicates that their flux 
is usually difficult to distinguish from that of the un- 
derlying basal flux. That is, the two components have 
about the same high temperature of 10-12 keV. 

Apparently random variations of several hours and 
lO's of percent amplitude occur often in the 7 Cas 
X-ray light curve. One long simultaneous observation 
with the Goddard High Resolution Spectrograph (GHRS) 
on the Hubble Space Telescope established the corre- 
l ation of small-amplitude UV con ti nuum fluctuations 



( Smith. Robinson, fc Corbet) 119981 : ISmith fc Robinson 
19991 ) and the concomitant absorptions in either lines aris- 
ing from either lower or higher than expected ionization 
and excitation states (SR99). 

Long-term variations have also been reported with cy- 
cle lengths of 50-91 days and amplitudes of a factor of 
three (Robinson, Smith, fc Henry 2002, "RSH02"). These 
are so far well correlated with corresponding cycles in op- 
tical bands (Johnson B and V) . Although the optical am- 
plitudes are 100 times smaller than the X-ray ones they 
are still much larger than in terms of luminosity output 
and therefore cannot be the result of reprocessing of X- 
rays. Rather, a common mechanism appears to mediate 
both variabilities. 

Modulations with several hour s timescales we r e sug - 



gested by previous X-ray satellites. iFrontera et al.l (|l987r ) 
reported a modulation with timescale of ~ 1.67 hr 
from EXOSA T obse rvation on 1984 December 7, while 
Parmar et al.l (|l993f ) argued from a rcanalysis of this 
data that such an oscillation arise from statistical fiuc- 
tuations in the red noise spectrum of the source. Also, 
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Parmar et alj (jl993[ ) have not found periodic oscillation 
in other EXOSAT data obtained on 1984 December 25- 
26. A 



2.3 hr oscillation was suspected by|f 




from R OSAT observation on 1993 July 16-17 

bund weak evidence for a ^ 0.6 h and 2.6 hr in 
BeppoSAX observation of 7 Cas carried out on 1998 July 
20-21. None of these suspected detections could be found 
by subsequent monitorings of the star with the RXTE 
(e.g., Smith, Robinson, & Henry 2000). 

Pre-high resolution spectroscopy of 7 Cas has demon- 
strated that its X-ray emission is domi nated by a ther- 
mal c omponent with fcT ^ 10-12 keV (Murakami et al. 



1986 : IParmar et al.l[l99l iKubo et allll998t lOwens et al 



19991) . These spectra exhibited FeXXV and FeXXVILya 
K line strengths consistent with this temperature, but 
they are consistent only with a subsolar Fe abundance 
of 0.2-0.4 Z0. 

Prior to this paper, the only high dispersion spec- 
trum of 7 Cas was obtained on 2001 August 10 by the 
Chandra High Energy Trans mission Grating (HTEG) dis- 
cussed bv lSmith etall (|2004l) (hereafter "S04"). This spec- 
trum showed a complex structure caused by plasma ra- 
diating in at least 3-4 temperatures ranging from about 
12 keV to about 0.15 keV. The hot component in turn con- 
sisted of two subcomponents with very different (~10^^ 
cm~^ and 3x10^^ cm~^) column densities. While the FeK 
lines again indicated abundances of only O.22±O.O5Z0 
other Fe lines arising from the L-shell give abundances 
which were consistent with the solar value. The strengths 
of other lines were also consistent with solar abundances. 



3. Observations and data reduction 

7 Cas was observed by the XMM-Newton X-ray observa- 
tory on 2004 February 5 during about 68 ks in the rev- 
olution 762 (Obsid 0201220101). This observation was 
performed with the EPIC pn camera running in the fast 
timing mode (timing resolution of 0.03 ms), connected to 
the thick optical filter, and with the high spectral res- 
olution RGSl and RGS2 cameras. The optical monitor 
was blocked as usual for bright optical sources, and the 
central CCDs in the MOSl/2 cameras were switched off 
in order to avoid overloading the telemetry. Data reduc- 
tion has been made using the Science Analysis System 
(SAS) software vS.O.l. All data were reprocessed using the 
pipeline epproc (for EPIC pn camera) and RGSPROC (for 
RGS camera s) tasks. For the tim ing analysis, we use the 
Rayleigh (Buccheri et al.lll983l ) and the Scargle/Midas 
( Scarglelll98a^ peridograms, and the Xrono^ package. For 
spectral fits, the XSPECH software vll.3.2 was applied. 
The EPIC pn data obtained in the fast timing mode have 
a doubtful calibration and increased noise especially at 
softer energies (< 0.5 keV: iGuainazzi 2008 ') and, to be on 
the safe side, timing analysis was restricted to broadbands 
at the 0.8-10 keV energy range without relevant loss of in- 



formations. For spectral studies, since we have high qual- 
ity spectra from RGSl and RGS2 cameras covering the 
soft X-rays (6.2-38A; or ~ 0.3-2 keV), we use the EPIC 
pn data as complementary data covering the hard part of 
the spectrum (1.25-4.lA; or ~ 3-lOkeV). We noticed that 
the inclusion of the low energy part of the EPIC pn data 
produces systematic residuals that can be connected to 
calibration uncertainties. Because of the high count rate 
of 7 Cas, the inclusion of times with slight background 
flares during the XMM-Newton observation has little im- 
pact on its EPIC pn and RGS light curves. However, we 
opted for excluding these times in the spectral analysis. 
The resulting exposure times in the flare-free good time 
intervals were 48 ks for the RGSl and RGS2 cameras and 
51.2 ks for EPIC pn. 

We reinvestigated the Chandra data of 7 Cas ob- 
tained on 2001 August 10 with the HETG during about 
52 ks (ObsI D 1895 ) with the CIAOv4.0, reported by 



^ http:/ /heasarc. nasa.gov/docs/xanadu/xronos/xronos. html 
^ http:/ /heasarc. nasa.gov/docs/xanadu/xspec/index. html 



Smith et al. ()2004l) . Our purpose was to compare the 
Chandra and XMM-Newton spectra of 7 Cas from the 
same analysis techniques and thus to insure that any dif- 
ferences between the derived parameters are not due to 
the model fitting programs. Our reanalysis of the Chandra 
spectrum with XSPEC and families of mekal models gives 
results consistent with those reported by S04, obtained 
from Sherpa and APEC models. 

4. Timing studies 

The soft (0.8-2 keV) and hard (2-10 keV) light curves of 
7 Cas are marked for strong variabilities on timescales 
ranging from a few to thousand of seconds (Fig. [T]). On 
long timescales the variation attains ~80% of the mean 
flux. Superposed on this global modulation are ubiqui- 
tous flare-like events. We note in particular that the soft 
X-ray light curve responds faithfully to the rapid flares 
of the high energy light curve. In addition, the hardness 
of the source (Fig. [TJ;) usually exhibits nearly the same 
variations as the integrated flux. For example, the formal 
slope of the hardness-total count rate curve from our ob- 
servations (not shown) is only (5.2±2.9) x 10^**. 

These findings are similar to those reported by SRC98, 
though RSOO found that shot fluxes tended to exhibit a 
lower (softer) ratio. 

Figure [5]a, exhibits the power spectrum of the EPIC pn 
events from the combined spectral bands shown in Fig.[TJ 
A key point from this figure is that at frequencies above 
about 0.003 Hz the slope of the distribution is significantly 
steeper than -1. This was also noted by SRC98 and RSOO, 
who found slopes of -1.23 and -1.36 in 1996 and 1998, 
respectively. The difference between those values was al- 
ready marginally statistically significant, and the current 
result is certainly even more significantly different from 
the 1996 result. As noted by RSOO, the difference between 
these slopes and -1 is most likely due to a relative preva- 
lence of strong, longer-lived flares as compared to rapid 
short ones in these data. RSOO reported 1// slopes in their 
1998 light curves that were intermediate between the 1996 
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Fig. 1. RGS1-I-RGS2 light curve combining 1st and 2nd orders and the RGSl background light curve (a), and EPIC 
pn light curves in the 0.8-2 keV (b) and 2-10 keV (c) X-ray bands, and the respective hardness ratio (d). Time bins 
of 100 s. Upper limits for error bars at 1 ct are plotted in the left corner of the figures. The horizontal lines in panel 
(b) represent the thresholds adopted to define the high (H) and low (L) state spectra (see Section [5?5|) . 
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Fig. 2. Power spectrum. Top: from EPIC pn events at 
0.8-10 keV, rebinned as a geometrical series. Bottom: from 
0.8-10 keV EPIC pn data binned to 100 s; the dashed lines 
represent the confidence levels. 



slope of -1 and the present steeper one. The explanation 
for the knee at about 0.003 Hz is again likely to be caused 
by the absence of discrete flares with timescales longer 
than a few minutes. 
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Fig. 3. Autocorrelation from 0.8-10 keV light curve 
binned to 500 s. PI refers to result of Fig. Oj- 



The break is due to the dominance of apparently ran- 
dom variations that often occur in the X-ray light curve on 
timescales of about a half hour or longer. At least some 
of these appear to b e connected with v ariations in the 
UV continuum, which Smith. Robinson, fc Corbg^ ()l998f ) 
have associated with the partial occultation of the star by 
rotationally adverted translucent corotating clouds. 

Another interesting feature emerging from the auto- 
correlation analysis of a light curve obtained by the RXTE 
satellite (Fig.[3|) was the reap pearance of period i c flux 
lulls. Following the procedure of iRobinson fc Smith ( 2000l ) 
we have cross-correlated both the fluxes in our XMM- 
Newton light curve and also the reciprocals of these fluxes. 
The latter exercise was carried out because RSOO had 
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4. The high resolution fluxed spectrum of 7 Cas, from the RGS1-I-RGS2 XMM-Newton cameras combining 1st 
2nd orders, and from the HEG+MEG Chandra cameras combining all four first-order spectra. 



found that the autocorrelation of the direct fluxes pro- 
duced no distinguishable features, whereas the autocor- 
relation of the reciprocal fluxes produced dramatic vari- 
ations with lags of spacing of 7.5 hours and higher mul- 
tiples. These authors reasoned that this different behav- 
ior could come about if X-ray active centers at various 
latitudes of the rotating star would not produce marked 
features in the autocorrelation curve. However, features 
would be observed from the curve generated from the re- 
ciprocal fluxes if only a few longitudes of the star were 
not covered with X-ray activity. Our results, shown in 
Fig. [3] for the direct flux, produced rather similar curves 
for both the direct and reciprocal fluxes, indicating that 
the distribution of X-ray centers had rearranged them- 
selves dramatically between the two epochs of observa- 
tion. The curve itself shows evenly spaced maxima every 
10156 s (2.8 hours; Fig.j^b)- On a subsequent study of six 
RXTE light curves of 7 Cas, inson. Smith, fc Henrvl 
(|2002[l discovered that autocorrelation peaks were most 
apparent in light curves associated with flare fluxes. These 
authors found that the features are centered at differ- 
ent time lags at different epochs, and occasionally did 
not occur at all. At various times these authors found 
spacings at 7.5 hours, 3.5 hours, and 5.8 hours. The 
2.8 hours spacing we now find appear to be a recur- 
rence of the same cycle length these authors found dur- 
ing a long RXTE observation in 2000 December. As 
Robinson. Smith, fc Henrv ( 20021 ) found, the autocorrela- 
tion peaks seem to be caused by an absence of flux (lulls) 
at different times. As such, they seem to hint at the exis- 
tence of X-ray relaxation cycles in the 7 Cas environment. 



In contrast. iLopes de Oliveira et al.l (|2007| ) found no such 
lulls in the light curve of HD 110432. 



5. Spectroscopic properties 

5.1. Comparison with previous Chandra HETGS 
spectrum 

We have obtained a XMM-Newton spectrum in order to 
investigate time-dependent differences and also to exploit 
the higher effective aperture of the system, albeit with 
lower resolution at higher energies. The XMM-Newton 
provides coverage over longer wavelengths (~ l-BSA) than 
an Chandra spectrum with the high and medium energy 
gratings ('- 1.5-25A). 

Our XMM-Newton spectrum shows the same general 
properties as the 2001 Chandra spectrum (Fig. Both 
spectra show continua consistent with a multi-component 
thermal model that includes a number of Lyman a lines 
of hydrogen-like ions, helium-like ions, as well as a few 
Lyman /3 lines. Fluorescent lines from lower ion stages of 
Fe and Si are also visible. The high energies are dominated 
by the presence of the "FeK complex" (FeXXVILya, 
FeXXV, and fluorescent FeK lines) in the range 1.7- 
2 A. The two spectra of this aggregate look similar, ex- 
cept that the fluorescence feature is weaker in the XMM- 
Newton spectrum. However, our RGS spectrum reveals 
many more details than the Chandra one because of the 
XMM-Newton^s larger effective aperture as well as an ex- 
tension to longer wavelengths and because in 2004 the 
soft X-rays of 7 Cas were attenuated far less than in 2001. 
These circumstances allowed us to probe to lower tern- 
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peratures and to better judge the number of components 
required to fit both hues and continua in the soft X-ray 
region. 

5.2. Selection of multi-component modeling 

Our models are the results of fitting both the continua and 
the line strengths in a multivariable solution with either 
the mekal formulation for optically thin plasmas, or a re- 
lated model, vmekal, in which individual elemental abun- 
dances can be determined independently. A Gaussian line 
was included to describe the fluorescence feature arising 
from low-ion stages of Fe at 1.94 A (6.4 keV), not incorpo- 
rated in the mekal code. In all cases we adopted the phabs 
code to describe the photoelectric absorption. When re- 
quired, we_assumed_the Hipparcos distance of 188 pc for 
7 Cas ( PerrvmanlllQQTl) . 

The high quality of our EPIC pn and RGS spectra al- 
lowed us to determine in detail the thermal components 
required to fit the line and continuum fluxes. This was 
especially important in allowing us to test whether the 
plasma has a continuous range of temperatures. Initial 
explorations showed that contributions from 3 or 4 tem- 
peratures are necessary. In addition, otherwise mediocre 
fits to the continua (especially the soft continuum) re- 
quired us to consider a two-absorption column model. We 
found first the same hot temperature that many other 
investigators have reported. We will refer to this compo- 
nent as kTq. This component has a 12keV temperature 
and dominates the total X-ray flux. The carbon and ni- 
trogen line spectra are emitted by a cool plasma with a 
temperature fcTi of ~ 0.1 keV. The unabsorbed luminosity 
of this component, integrated over the range 0.2-12 keV, 
is 6.8-8 X 10'^^ erg s~^. One or two intermediate (warm) 
temperatures are likewise required to explain the pres- 
ence of lines of oxygen, neon, magnesium, silicon, a few 
Fe ions with L-shell configurations (Fe XVII- XXIV), and 
the Fe XXV line (Figures H and [H) . 

Before refining the initial 3 and 4 component models 
further, we attempted to determine whether these multi- 
ple thermal components were discrete or could be part of a 
continuous distribution, as might be expected in a single, 
thermally differentiated plasma. To continue tests using 
this initial set of simple models, we froze all but one tem- 
perature in our 4-T model and computed models for the 
remaining temperature specified over graduatedsteps. We 
made movies of these results overplotted with the obser- 
vations and determined those values of the scanned tem- 
perature which allowed us to judge the best agreement 
with the observations. This technique also allowed us to 
note whether predicted lines are not observed as well as to 
test whether the principal plasma components are essen- 
tially monothermal. FigurelHl demonstrates this for fcT2 in 
M2 (4-T), which is necessarily confined within the range 
0.5 and IkeV. This plot shows the absence of Fe XXIII 
and strong FeXXIV lines, as marked on our IO-I3A seg- 
ment of our RGS2 spectrum (RGSl shows a gap at this 



high state 
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Fig. 5. The FeKa complex and its best fit determined 
from the model bremsstrahlung -|- 3 Gaussian lines for the 
low and high-state spectra (bottom and top, respectively) 
and the integrated spectrum. See details in Section [5.51 



wavelength). These lines would be stronger visible if the 
temperature were in the range 0.7-0.9 keV, and yet they 
are not seen. This fact eliminates the possibility that fcT2 
is a distribution of temperatures extending to higher tem- 
peratures than 0.7 keV. On the lower bound our models 
predict the presence of the Fe XVIII 14.2 A line for tem- 
peratures <0.5keV, and this is not visible either. In ad- 
dition, the abrupt disappearance of components from the 
OVII complex in the models shows that the absence of 
intermediate temperature plasma extends from 0.5 keV to 
about 0.2 keV. Altogether these diagnostics indicate that 
the fcT2 component in the 7 Cas environment is limited to 
a very narrow range of temperatures centered near 0.6 keV 
and that it is distinct from a lower or higher plasma com- 
ponent. This statement can probably be extended to the 
spatial separation of the warm and cool plasmas as well. 

Proceeding with this analysis, the weak MgXII Lya 
line hints the existence of a warm kT^ component at ^ 
2keV, although some of our models have managed to 
achieve fits without it. The absence of yet other Fe L lines 
indicates that there is likely no plasma emission in the 
range of 2.5 to at least 5keV. Also, the presence of the 
FeXXVI Lya feature, the main line diagnostic in the hard 
X-ray spectrum of 7 Cas, requires a temperature of at least 
8 keV. The Fe XXV Lya line requires a value significantly 
lower than this limit, and this means that it cannot be 
formed in the fcTg plasma component alone. In FigurelS] 
a feature due to NeX Lyman/? is present. This feature, 
visible in both the individual RGSl and RGS2 spectra, 
cannot be easily fit with our thermal equilibrium {mekal) 
models. 

It is quite possible that the kTq component consists 
of some distribution of temperatures around a mean value 
of 12-13 keV. This inference is supported by the occa- 
sional color variations in the high energy light curve of 
7 Cas (SRC98), suggesting the presence of hot many sep- 
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- kT2 = 0,40 keV 
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Channel vvtiVelengLh (A) 

Fig. 6. The observed RGS2 spectrum in the lO-lSA re- 
gion for three trial values of kT2 bracketing our model M2, 
in each case offset vertically for convenience. This figure 
depicts the Ne X Lya and predicted positions of Fe lines 
arising from several ions, and the preferred temperature 
kT = 0.63 keV for fcT2. These predictions allow us to place 
limits on the lower limit to the range of the warm compo- 
nent, kT2 (see text). The emission bump at 10.23 A arises 
from NeX Ly /?, and it is not fit in the range of tempera- 
tures shown. 

arated exploding gas volumes having a small range of tem- 
peratures. 

The upshot of these considerations the presence and 
absence of various lines in this spectrum forces the con- 
clusion that the plasma does not exhibit a continuous 
Differential Emission Measure (DEM). Rather, there are 
gaps over the ranges of about 2.5-8 keV, 0.7-1 keV, and 
0.2-0.5 keV. S04 had suspected that a continuous DEM 
was inappropriate but could not state it conclusively be- 
cause of the pronounced photoelectric absorption of the 
Chandra soft-energy spectrum. 

5.3. Model fitting 

The above procedures allowed us to refine our initial 3-T 
and 4-T models, utilizing both line and continuum fiux 
information. Further attempts to search for a continuous 
Differential Emission Measure from the continuum alone, 
such as with a cooling flow cemekl led to significantly de- 
graded fits even in the continuumlfl 

^ The cemekl/ cevmkl code describes a multi-temperature 
plasma based on the mekal code, in which the emission mea- 
sures of the plasma scale with their temperature as (T/Tmax)"; 
a = 1 corresponds to the adiabatic case. 



The fits converged for fcT^ax ~ 42keV and xl = 2.4. 
Freezing fcTmax to 12keV produced an unacceptably high 

= 6.4. We next tried a composite model cemekl + 
mekal. This resulted in temperatures of fcT^ax ~ 12 keV 
and fcT ~ 0.9 keV for each component, respectively, and 
a xt = 1.75. This model failed to describe the contin- 
uum at high energies and underpredicted the following 
features: the FeXXV line in the FeKa complex, the Fe L- 
shell lines, OV III Ly^, the fir triplet of O VII, N VII Lya, 
and NeX Lya. Allowing the cemekal a parameter to vary 
did not improve the fits. 

Table [T] shows the results of our analysis in terms of 
three basic models. Columns 2, 3, and 4 list the values de- 
termined for the temperature, column density, the mod- 
eled unattenuated fiux, the resulting percentage of the to- 
tal X-ray flux, and emission measures for the three and 
four component models, and the Fe abundance derived 
from the FeXXV and Fe XXVI Lya ions. In all cases the 
high temperature component we designate as fcTg was 
found to lie near 12-14 keV. Column 2 of Table[T]shows the 
solution for a 3-T model, "(Ml)," that includes a cool and 
warm components (fcTi and fcT2, respectively) as well as 
the dominant one, kTq. The cool/warm component fluxes 
were attenuated by a column Nh^ while the hot one was 
attenuated by Nh^ ■ Our 4-T model (column 3 of TableHJ 
"M2" ) gives the solution for the addition of a warm com- 
ponent having fcTa and also attenuated by Nh^^ • The final 
column of the table, "(M3)," is the 3-T solution again. 
However, this time the hot component was attenuated by 
both Nfjj^ and Nhi, columns. In this model, the dominant 
hot subcomponent (fcTg), affected by the individual Nfjb 
column, contributes to ^ 76% of the total flux. The second 
hot subcomponent (fcT'g), warm (fcT2), and cool (fcTi) 
components were affected by a common Nh^ absorption 
column. Figures [7] and m show the unfolded spectrum for 
each model, allowing the reader to judge which of the var- 
ious observed lines is formed in which plasma component. 
In our solutions we found that Nh.^ ~ 2x10^^ cm^^ and 
~ 2-3x10^" cm~^. By comparison, the UV and Ha 
determined ISM column density to 7 Cas from the litera- 
ture is a scant 1-2x10^° cm~^, which leaves only a little 
room for absorption of the soft X-rays within the source. 

The 4-T fit we found included components fcTi, kT2, 
fcT3, and kTq having values of about 14keV, 2.4keV, 
0.6 keV, and 0.1 keV, respectively. For the 3-T models, the 
fcT3 component was omitted, but this omission must be 
compensated for in the Fe XXVI Lya strength by decreas- 
ing the temperature from 14 to 12.5 kcV. 

It is worth noting that any of our models give us the 
desired 1- The residuals are mostly in the lines, not 

in the continuum, and the relatively high - although even 
acceptable - values of x^ (Table [IJ are due in part to an 
inadequate description of the individual line profiles. 

A key point to take away from the Table [T] is that al- 
though 84-90% of the total fiux is radiated by the kTq 
(or kTq+kTg) component(s), almost all the lines except 
for the FeK complex is contributed by the cooler A;Ti_3. 
The hot kTq component therefore provides no kinematic 
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Table 1. Spectral fit results using the 3-T and 4-T models. 





hot component, 1-col. 


hot component, 2-col. 




0- 1 


A T 


3-T 




(Ml) 


(M2) 


(M3) 


Nh, (10^2 cm-2) 


— 0.03 




0.24t°:°* 


fcTi (iceV) 
/ti (erg cm~^ s"""^) 


0-10-0.01 


n 1 1 +0.01 


O.lllo.ol 


i.DXiu yl.o/0) 


1 Q\^in~^^ f^A^A 
~ i.yxiu yo.^/o) 


~ 1.8x10"" (8.0%) 


EMti (10^^ cm-3) 


~ 0.5 (11.5%) 


~ 0.7 (15.4%) 


~ 0.6 (12.9%) 


kT2 (keV) 
/t2 (erg cm"^ s~^) 


0-t>4_o 04 


n RA+0 03 


0.62t°:°^ 


~ 4.UXiU t^i.o/oj 


~ 4.zXiU t^i.y/oj 


~ 4.0x10-12 (1.8%) 


EMt2 (10^^ cm-=^) 


~ 0.05 (1.1%) 


~ 0.06 (1.3%) 


~ 0.06 (1.3%) 


fcTs (keV) 
/t3 (erg cm~2 s~^) 




2.40+^-2^ 

— O.ZD 






~ 1.2x10"" (5.3%) 




EMt3 (10^^ cm-3) 




~ 0.2 (4.4%) 




kTQ, (keV) 
/tq, (erg s"^) 






12.38«-j° (*' 






~ 3.2x10-" (14.3%) 


EMtq, (10^^ cm-^') 






~ 0.6 (12.9%) 


Nh^ (W cm-2) 




022+° °°? 

^•^^'^ — 0.002 




fcTg (keV) 
/tq (erg cm"^ s"^) 


19 cc+O.Sl 
lZ.OO_o 29 


IzL Q9+0-59 


12.38 <*' 


~ 2.0x10"^° (90.9%) 


- 1.9xl0~^° (84.4%) 


~ 1.7x10-1" (75.9%) 


EMt<j (10^^ cm"^) 


~ 3.8 (87.4%) 


~ 3.6 (78.9%) 


~ 3.4 (72.9%) 






>-'-ro_Q Qj^ 


]0+° °2 


l2+° °i 


Line (keV) 


6.4'' 


6.4'' 


6.4'' 


O-Line (keV) 


0.01 


0.01 


U.U»_Q 03 


ftot (erg cm"^ s"^) 


~ 2.2x10"^° 


~ 2.2x10-1° 


~ 2.3x10-1° 


EMtot (10^5 cm-^') 


~ 4.4 


~ 4.6 


~ 4.7 


x'/d.o.f.^ 


1.47/1532 


1.44/1529 


1.46/1530 



^ degrees of freedom; ^ frozen parameter. 

Notes: Fluxes are given unabsorbed in the 0.2-12 keV energy band. In parentesis are the fluxes and EM in % of the total 
values. Plasmas with solar abundances, except for the [Fe] of the hot component. Quoted errors are at the 90% confidence 
level. Base model: Nh,*(Ti-|-T2-|-T3-|-Tq/)-|-Nhi,*(Tq); see Section [53] for details. 



information. However, the possibility of multiple absorp- 
tion columns and the presence of the fluorescence features 
of Fe and Si do provide potential geometrical descriptions 
for the emitting volumes. 

Although S04's solution for the 2001 spectrum was 
similar in that it required a 4 thermal component fit, there 
were some small differences in detail that are significantly 
different. One of the clearest differences is with the pre- 
cise temperature of a warm component (their "fcTa" 
0.4 keV cf. the kT2 ^ 0.6 keV herein). A small increase 
in the OVIII/OVII Lyman a line ratio for the XMM- 
Newton spectrum is indicative of the higher value we find. 
If a warm temperature is a quasi-permanent feature of the 
X-ray spectrum, it has definitely shifted from 0.4 keV to 
0.6 keV between 2001 and 2004. The hot (- 12-14 keV) 
and cool (~ 0.1 keV) component temperatures overlap. 
Although we cannot be sure if ~ 2keV plasma existed 
in 2004, its possible presence is consistent with emission 
powered by a similar thermal component in the 2001 spec- 
trum in their 2001 Chandra observation. 



As already noted, one of most important differences 
from earlier results is that no high absorption column was 
needed to fit this spectrum. S04 found that the Chandra 
spectrum required a fitting with two hot plasma absorp- 
tion columns. The first contributed 10-30% of the 12 keV 
component and was attenuated by an high absorption col- 
umn of 10^^ cm"^. The second component, which is the 
dominant fiux contributor to the observed spectrum at 
high energies, was attenuated by a column of ~ SxlO^'^ 
cm-^. In contrast, most of the X-ray flux observed by 
XMM-Newton is affected by one dominant absorption col- 
umn equivalent to Nh ~ 2-3x10^° cm-^, as derived from 
our 3-T and 4-T model with a single absorption column. 
However, for our models we were only able to fit the con- 
tinua of our combined EPIC pn and RGS spectra with two 
absorption columns. Although the absorption of the Nfjb 
column, at ~ 2-3x10^*^ cm-^, is very low, it has much in- 
fiuence because it affects most (62-96%) of the total emis- 
sion (see Table [1]). The other absorption column has only 
a minor influence in absorbing the flux even though it is 
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Fig. 7. Unfolded spectrum from RGSl and RGS2 for each model in Table [T] The blue, green, cyan (light blue) and 
red lines correspond to the fcTi, fcT2, fcTa and kTg (and kTgi), respectively, while the black lines correspond to the 
composite model. 



higher by a factor of 7-10 (Nh^ ^ 2x10^^ cm^^). We note 
that the temperatures derived for each component in the 
models with one, as opposed to two, absorption columns 
are consistent with one another. The two column model 
also improved the descriptions of the lines formed at low 
energies, such as NeXLya, and OVIIILya. The improve- 
ment of the fit with the second column over a model with a 
single column is supported by the F-test statistic, stating 
with a probability of only 4.7x10"^ that the derived de- 
crease in xt associated with the second column column is 
due to chance. The important conclusion from this is that 
the absence of a high absorption column as derived from 
the XMM- Newton observation suggests a dramatic rear- 
rangement of cold circumstellar gas in the vicinity of the 
X-ray sources between 2001 August and 2004 February. 
Contrary to what is arguable for the 2001 epoch, it can- 
not be said that in 2004 the circumstellar disk was in front 
of a significant fraction of the soft X-ray source(s). 

In order to estimate the origin of the long-wavelength 
lines in our models, we constructed a lineless contin- 
uum by excluding of the spectral regions lines in the 10- 
I7A range and fit the resulting pure continuum spec- 
trum to a bremsstrahlung model. This model gave a 



value a slightly hotter temperature 13.2 keV for our hot 
component with a xt of 1.17. The addition of a second 
bremsstrahlung component did not improve the fits (/cT 
0.04 keV; xl of 1.17). This result suggests that the de- 
termination of the low component temperatures, fcTi and 
kT2, is strongly driven by the lines alone. 

Although all of our discussion has taken place in the 
context of equilibrium models, we agree with Smith et al. 
(2004), who suggested that the warm and/or cool plasma 
components could also result from nonequilibrium ioniza- 
tion processes, e.g. from the sudden impact of gas parcels 
ejected by flares into a dense stationary medium such as 
the Be disk. We suspect this could also be the reason why 
we observe Lyman /3 transitions, such as the Ly/? tran- 
sitions of NeX and OVIII, that are not easily fit with 
equilibrium mekal models. A similar explanation might be 
also important in understanding the Fe XXVI Ly/3 line at 
8.3 keV of HD 110432, reported bv lLooes de Oliveira et"ar 



(|2007l ). and which could not be fit with standard models 
for the ranges of temperature considered for this feature. 
The strong Lyman /3 features hint at excitation tempera- 
tures in these plasmas that are higher than the ionization 
temperatures we found in Table [TJ 
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Fig. 8. Unfolded spectrum from EPIC pn, centered at the 
FeKa complex for each model in Table [T] (see Fig. [7] for 
details). The dotted lines refer to the Gaussian line at 
1.94A(6.4keV). 



5.4. Abundances and anomalous line strengths 

For the most part the line strengths are consistent with 
solar abundances in our 3-T and 4-T models. The first 
and most conspicuous exception to this general result is 
the finding from Table[T] that the Fe abundance from K- 
electron Fe ions is 0.12±0.02Zq. This is a factor of two 
lower than any other studies in the literature (the S04 re- 
sult was 0.24±0.02 Zq). Because K-shell ion abundances 
of 0.1 Zq and 0.24 Z© are significantly different from one 
another, our result clarifies in a new way that the Fe abun- 
dance derived from K-shell ions changes with time. Thus, 
the K-shell result may not reflect a global elemental abun- 
dance. 

As also found by S04, the Fe abundance from L-shell 
ions is significantly higher than the Fe abundance from 
K-shell ions. We pursued the investigation of the Fe abun- 
dance from the Fe L-shell lines with the models Ml, M2 
and M3. In a new set of models we fixed other variables 
and solved for the Fe abundance of the soft temperature 
components fcTi and kT2. Our models suggest that the 
Fe L-shell abundance could be 0.4 to 1.3xZq. This exer- 
cise verifies the distinctness of the Fe abundance arising 
from the L-shell lines on one hand and the K-shell lines on 
the other. For completeness note that a similar K- and L- 
shell F e anomaly was reported bv lLopes de Oliveira et al. 
(|2007^ for HD 110432. 

The second apparent departure from solar abun- 
dances comes from the enhanced strengths of the lines 
of hydrogen and helium-like ions of neon and nitrogen. 
The formal abundances we find from our models are 



Zn=3.961[]J^xZn,0 and ZNe=2.63l:J]:^^xZNc,0, assum- 
ing solar abundance for the other metals in the softest 
plasma and allowing for line broadening and bulk ve- 
locities. However, the allowed range of Ne abundances 
depends on the number and nature of the parameters 
left free in the fit, the most important being the as- 
sumed Fe L-shell abundance and the line broadening ve- 
locity. The hottest component dominates the NeX line 
flux with a smaller (^ 1/3) contribution from the T2 
thermal component. For instance, using bvapec instead 
of mekal and leaving Fe L-shell abundance free yields 
slightly lower abundances with ZNo=l-65lg'2gXZNo,0 and 
ZFeL-shoii=O.4i;to'iiXZFe,0- This cffect is due to the pres- 
ence of a FeL-shell forest around the NeX La line and to 
the fact that the continuum level depends on the descrip- 
tion of the FeL-shell lines. Although the limits for the 
abundances are not well determined, there is a clear ex- 
cess above the solar values for N and Ne. The NeX line 
peak velocity is in the range of -110 to -1-230 kms~^ and 
displays a broadening velocity comparable to that of other 
bright emission lines. 

A similar analysis was applied to the Chandra spec- 
trum in order to determine whether the Ne enhancement 
was present or not in 2001 - it is not possible to expand 
the investigation to N abundance due to the low SNR 
of the NVIILa line. The non- velocity broadened mekal 
model yields Znb = O.75±O.19xZNe.0 consistent with the 
value obtained by Smith et al.l ( 20041 ) . However, since the 
observed NeX line apppears significantly wider than mod- 
eled, we used the bvapec model to take into account such 
a broadening. The larger NeX EW impacts the determi- 
nation of the abundance yielding ZNc=l-52l^g'24xZNe,0 
(or Znc = l-2Olg'24xZNo.0 ff the oxygen and iron L-shell 
abundances are left free, in which case we obtain ZFoL-shcii 
— O.46j;g QgXZFe,0)- The broadened line profile fits the ob- 
served Ne line much better, with peak velocities then in 
the range of -60 to 260 kms~^ and broadening velocities 
of ^ 500 ± 200 kms^^. These values are consistent with 
those derived from the RGS spectrum. 

We note carefully that whereas the fluxes in these lines 
scale linearly with abundance, the equivalent widths, as 
formally defined with respect to the neighboring contin- 
uum fiux, scale far more slowly because of the contribution 
to the local bound-free opacities from ions in a N-Ne rich 
plasma. Thus, we found an increase of only 35% with re- 
spect to the equivalent width measured in the Chandra 
spectrum. Our models with XSPEC validate this mild in- 
crease in equivalent width with abundance in detail. Given 
this reality, we anticipate that the true abundance errors 
are larger than those XSPEC computes based on photon 
statistics. Nonetheless, the anomalous excess cannot be 
discounted. Moreover, the possibility that the X-ray envi- 
ronment could be nitrogen-rich by a factor of 3-4 is unre- 
markable because enhancements are already a hallmark of 
massive stars evolving off the main sequence, although the 



reasons for this are still under discussion (jHunter et al 
l2008i) . However, the apparent neon enhancement is of 
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much greater interest because in a stellar evolution con- 
text neon is produced by carbon burning shortly before 
supernova detonation or in the interiors of some white 
dwarfs, and dredging to the surface of such elements is 
possible by various mixing processes. 

Suggested alternatives to the straightforward abun- 
dance interpretation of the line strengths are the following: 
(i) instrumental artifact or a cosmic ray; (ii) an inappro- 
priate temperature used in the modeling; and (iii) the line 
is strengthened by a microturbulent like broadening. Each 
of these possibilities may be dismissed in turn, (i) is un- 
likely because separate spectrum extractions from the pos- 
itive and negative detector halves show consistent profiles. 
Likewise, we may rule out (ii) because a temperature of 
0.6 keV is already ideal for NeX formation, and a stronger 
feature cannot be produced given a standard abundance. 
Possiblity (iii) would require a large optical depth and 
strengthening of the line. However, even though the col- 
umn density is probably high enough for photons to expe- 
rience more than one mean free path through their transits 
across the medium, they are nonetheless scattered coher- 
ently. Therefore the line widths are unaltered by these 
histories. With these possibilities ruled out, we are forced 
to conclude that the neon and nitrogen abundances in the 
X-ray plasmas of 7 Cas are high. 

In considering the anomalous Fe abundance derived 
from Fe-K lines, but normal abundance from Fe-L lines in 
the 7 Cas spectrum, S04 suggested that this could arise 
from an inverse FIP (first ionization potential) effect in 
the 7 Cas environment that is similar to that found in 
the coronae of the Sun and other ma gnetically and X-ray 
active cool stars, such as AB Dor ([Gudel et al. 1 12003,1 . 
Although the details are still unclear, it appears that in 
a magnetic plasma low density environment differential 
ponderomotive forces, produced by wave heating and de- 
pendent on their first ionization potential, can prevent cer- 
tain ions from migrating across a plasma and resulting i n 



Table 2. Parameters of the emission lines of the FeKa 
complex from a bremsstrahlung -I- 3 Gaussian lines model. 



an altered measured chemical abundance ([Laming 1120041 ). 
Whether such a process is actually active in the environ- 
ment of 7 Cas, let alone whether it extends to Ne^+ ions, 
must be considered speculative. Its attractiveness lies in 
its potential to explain the different Fe abundances de- 
rived from K-shell and L-shell ions. Similarly, now in the 
framework of the accretion model, it is not clear, first, how 
Ne could be preferentially ejected from the white dwarf at- 
mosphere or, second, how Ne could be enhanced without 
enhancing at the same time C, O, and Mg abundances. 



5.5. The Fe Ka complex 

In evaluating the equivalent widths (EW) of each emis- 
sion iron line of the FeK complex, we used the 5-10 keV 
photons acquired during low background phases of the 
XMM-A''ewton satellite's orbit. Our approach was to limit 
consideration of spectra accumulated only during low (L) 
and high (H) flux states observed in 7 Cas (defined as the 
lowest and highest one-third fluxes of the total distribution 
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Total flux in line, in units of photons cm~^ s~^. Notes: 
Quoted errors are at 90% confidence level. The tempera- 
ture of the bremsstrahlung model at 1.2-4.1 A (3-10 keV) 
converges to fcT = 13.36(-|-0.21/-0.25) keV. 



in our light curve) and also during some 1000 time win- 
dows in which the detector background was low (see Fig.[T]- 
b) . Finally, we applied an absorbed bremsstrahlung model 
in order to describe the underlying continuum and three 
Gaussian lines to account the iron lines. Table [5] and Fig. 
[5] show the results. The measured centroid energy of each 
modeled Gaussian line is, within the errors, in agreement 
with the theoretical values for the fluorescent, helium-, 
and hydrogen-like components of the FeKa complex. It 
is possible that there may be a slight inverse sensitivity 
of the fluorescence strength (relative to the Fe XXVI line) , 
though this inference is of marginal statistical significance. 

In our EPIC pn spectrum the strength of the fluo- 
rescence feature (EW ^ -10 mA) is weaker than whe n 
Chandra observed it (EW - -19mA; ISmith et al.ll2004l ). 
in line with the fact that the attenuation of soft X-rays 
by cold matter is less. This suggests that the fluorescence 
emission feature is formed at least partially by the same 
medium that aborbs the soft-X ray flux. 

Curiously, the analysis in the low and high-state spec- 
tra reveals that the FeK fluorescence feature may be 
marginally stronger in the low flux case. Although such in- 
vestigation could not be expanded to a quantitative anal- 
ysis of all spectral range because the relatively limited 
signal-to-noise in the final spectra for RGSl/2, we sus- 
pect that there is weak evidence for a weakening of the 
SiK feature for the low flux spectrum. 



5. 6. Helium-like diagnostics of electron densities 

Our spectrum covers the regions of the He- like NelX, 
OVII and NVI complexes - but the NVI complex falls 
partially onto gaps in the CCD for RGSl and RGS2. 
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Each of these is comprised of a so-called fir (forbid- 
den/intercombination/resonance) line triplet (see Fig.|3]), 
and the ratio of their intensities can indicate whether the 
dominant excitation process for producing this triplet is 
coUisional or photoionization. For example, if collisions are 
dominant, the ratio G ^ {i + f) /r ^ 1 obtains, whereas 
if photoionizations dom i nate the exitation G ~ 4 (e.g. 
Porquet fc Dubaul I2OOOI : IPorquet et ai] l200l[ ). From the 
O VII fir complex in our XMM-Newton spectrum, we es- 
timate G ^ 0.9±0.1 for 7 Cas. The NVI rif components 
are poorly measured because they fall within the gaps of 
the RGS2 detector. Nonetheless, by computing their ratios 
we find nearly the same value, G ~ 0.7±0.3. This suggests 
that the plasma is in the classical domain of collisional 
dominance. As it happens, the fir ratio alone cannot dis- 
tinguish between quenching of the forbidden transition by 
collisions or by photoexcitations by a nearby strong UV 
source, such as the Be star. 



5. 7. Velocity broadening and shifts in warm plasma 
lines 

In Table [3] we list the observed centroid wavelength, the 
equivalent width and the flux of lines for those lines for 
which meaningful measurements could be made. The best 
observed lines in the soft X-ray spectrum, NeXLya and 
OVIIILyo; profiles are noticeably broadened. We evalu- 
ate the broadening of the lines by using the bapec model, 
a velocity- and thermally-broadened emission model from 
coUisionally-ionized diffuse gat|f|. The Gaussian sigma for 
velocity broadening converges to ~ 400 kms~^ for a model 
like Ml, M2 and M3 in TablelU replacing the mekalhy the 
bapec code. The line broadening was also measured on in- 
dividual lines, not only on the whole spectrum. For exam- 
ple, we found a velocity broadening of 300±70 kms~^ for 
the OVIII line. The RGS velocity is consistent with that 
reported by S04 (478±50 kms~^). We applied the same 
analysis to two RGS spectra of AB Dor acquired in 2000 
and 2006 (obsID 0126130201 and 0160363001 with up-to- 
date calibration) and found in each case that the broad- 
ening of the OVIII line was less than about 70 kms~^ at 
the 90% confidence level, proving that the lines of 7 Cas 
are indeed broadened. 

In order to check for velocity shift in the lines, we run 
the models in Table [T] leaving the redshift parameter free. 
We found a overall redshift of ~ 200 kms~^ - and values 
for the parameters and xl consistent with those shown 
in Table [T] - slightly larger than the typical error on the 
absolute wavelength scale of 7 mA, or ~ 140 kms~^. 

5.8. Do our cool or warm plasma component arise in a 
radiative wind? 

It has been believed for some time that X-ray emis- 
sion in O and early B stars arises in sites associated 
with shocks distributed in the stars radiative winds (e.g. 



Leutenegger et al.ll2006l:ICohenll2008l ). Analyses of HETG 
spectra of other stars in the same general region of the H-R 
Diagram, for example ofrSco (B0.2V), /3Cru (BO 5 III), 



and COp h (09 5 V) bv lCohenI (HqqI) , |Cohen 



and 



Waldron respectively, paint a complex picture. 



The Cohen et al. analyses show that the winds of at least 
some early B stars cannot be described by a standard pic- 
ture of wind-shocks, and indeed wind structures of at least 
T Sco and ( Oph may be influenced by the effects of chan- 
neling by magnetic fields in the regions surrounding the 
star. The temperatures of the gas are not yet known to 
vary but taken as a group range from about 3 to 25 MK. 
We note that the X-ray emission measures associated with 
winds in these stars have not been found to be variable so 
far. 

7 Cas has such a wind, as evidenced by the pro- 
files of doublets in UV resonance absorption lines of 
several ions, notably CIV, NV, and Si IV. These pro- 
files exhibit Discrete Absorption features at 1 000-1100 
kms~^ and edge velocites of -1800 kms~^ (IDoazan 



1982':'Ka per et al.lll 996': 'Sm ith. Robinson, Corbetlll998t 
iSmith. Robinson, fc Hatzes 19981 ). Our analvsis shows 
the presence of soft X-ray components whose tempera- 
tures, fluxes, and emission measures (Table [T]) are con- 
sistent with those of normal massive stars reported by 
(l2009h 



Walborn et al 



http://cxc.harvard.edu/atomdb/ 



for example. It is thus expected 
that at least part of such emission is emanating from the 
radiatively driven wind of 7 Cas. However, as can be seen 
in Fig. [3 the hottest component dominates also the soft 
part of the spectrum of 7 Cas and the contribution of a 
pure radiative wind is overshadowed by this hot compo- 
nent. Establishing the presence of a luminosity change in 
the soft components supporting a non-wind origin is a 
complicated issue because of the presence of the Be disk 
itself (which may shield a part of the wind) and because 
the characteristics of the 7 Cas' spectra in 2001 and 2004 
are different from each other. It is worth noting that the 
X-ray 0.2-12 keV (unabsorbed) fluxes from the soft (fcT 
■~ 0.1 keV) component in 2001 and 2004 are consistent 
with each other (< 1.7xl0~^^ ergs^^ from Chandra and 
1-2.7x10"" ergs"^ from XMM-Newton), while the flux 
from the fcT ^ 0.6 keV component has increased between 
2001 and 2004 (from 0.7-2.2 to 2.5-6x10-12 ergs"!). In 
conclusion, we cannot discard an additional contribution 
to the soft plasma in addition to that of a radiatively 
driven wind usually observed in massive stars. 



6. Conclusions 

We have reported the second high dispersion analysis of 
the X-ray spectrum of the X-ray anomalous B0.5e star 
7 Cas obtained by the XMM-Newton in 2004, and we find 
the following characteristics at this epoch: 
a) The emission is due to an optically thin, thermal 
medium comprised of 3 to 4 discrete components but dom- 
inated by a hot component having fcTg ~ 12-14 keV. The 
temperature of at least one component (fcT2 ~ 0.6 keV) 
has definitely shifted since the Chandra HTEG observa- 
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Ac 


EW 




Flux'' 




(A) 


(mA) 


(xlO-^) 


NeXLya'^ 


12.131[11] 


125[18] 


2.15[31] 


(1.26±0.16) 


OVIIILy/?'^ 


15.997[41] 


40[14] 


0.55[20] 


(0.31±0.16) 


OVIIILya'= 


18.9796[49] 


360 [28] 


3.93[31] 


(1.71±0.26) 


OVIIr'' 


2 1.60 [90] 


163[19] 


1.72 [20] 


(0.77±0.28) 


OVIH'' 


2 1.78 [92] 


149[19] 


1.58[20] 


(0.79±0.40) 


OYUf 


22.097^* 


<32 


<0.3 


(0.13±0.16) 


NVIILya= 


24.79[99] 


296[17] 


2.76[16] 


(0.61±0.03) 


NVIr'^ 


28.81 [90] 


<230 


<2 






29.09[95] 


<169 


<1.4 




NVI/= 


29.531'' 


<10 


< 0.08 




CVILya'' 


33. 736 [64] 


201 [29] 


1.37[20] 





^ Total flux in line, in units of photons cm ^s ^; in parenthesis are shown the values derived bv ISmith et all (|2004l ) from 
Chandra. ^ From RGSl. From RGS2. Frozen parameter. Quoted errors are at 90% confldence level. 



tion in 2001. It is possible that a kT^ 2.4keV component 
exists too. If so it may represent plasma with nearly the 
same temperature found by S04. Further, the presence of 
a cool component with fcTi ~ 0.1 keV is consistent with a 
value found by S04. This may or may not have the same 
cause as the wind-shocked fluxes emitted from other early 
B and Be stars. However, if so, it is a lower temperature 
than expected. Only part (at most) of the kT2 plasma 
could be produced in the wind, and if so the large line 
widths evidenced in this component indicate a larger tur- 
bulent broadening than is typical of winds in other early 
B stars. 

b) The hot kTg component appears stable and any dis- 
tribution of temperatures of individual sites around this 
mean value must be small. 

c) A subsolar abundance of iron is derived from the 
FeXXV and FeXXVILya features, in agreement with 
several other determinations. As S04 also found, this 
FeK abundance is significantly different of the abundance 
found from Fe-L ion lines. We have also discovered that 
[Fex] changes with time and was significantly higher in 
2004. 

d) The light curve of this star again shows ubiq- 
uitous, rapid flaring . Thi s was also found by 
iLopes de Oliveira et all (|2007t l in HD 110432. However, 
unlike HD 110432, color changes occurred only seldomly, 
and evidently not at all on the few hour timescale noted 
in HD 110432. 

e) Our light curve shows a quasi-periodic lull every 
2.8 hours in 2004, similar to the cyclical lulls of 3.5 hours, 
7-7.5 hours, and 5.8 hours noted by RSH02. These appear 
to occur in most epochs. 

f) A thick absorption column affected 25% of the hot 
component in August 2001, but it had disappeared by 
February 2004. 

g) Apparently variable FeK and perhaps SiK fluorescent 
features are present. These emissions correlate with an 
absorption column (point fj that attenuates soft X-rays. 



h) Broadening of warm (fcT2) component lines was notice- 
able at both epochs but may have increased marginally 
from 2001 (- 0.4 keV) to 2004 (- 0.6 keV). 

i) The strengths of lines of two ions each of N and Ne 
are underpredicted for XSPEC models with solar abun- 
dances. Upon consideration of alternative explanations, 
we have interpreted these as evidence of abundance en- 
hancements. However, we do not understand the cause of 
these enhancements. 

j) Given the identification of the Lyman /3 feature in a few 
ions formed in the warm and cool plasmas, there is more 
than a hint of nonequilibrium processes in the environ- 
ment of 7 Cas. This suggests formation in an (at most) 
intermediate density environment, which is clearly sepa- 
rate from the very high density plasma in which the flares 
are formed (Smith, Robinson, & Corbet 1998). 
k) The unabsorbed flux of 7 Cas at 0.2-12 keV in 2004 
from XMM-Newton 1.7-2.8x 10"^° ergs"^) is consis- 
tent with the value observed in 2001 from Chandra (~ 
2-3.1x10^1° ergs-i). 

In several of these respects the X-ray behavior of 7 Cas 
seems different from the one other analog studied in depth, 
HD 110432. In keeping with the differences of HD 110432, 
and indeed even the changing properties with time, we 
refer to this star as having "personality" - hence the title 
of this paper. 

Several changes we have noted are as remarkable as 
they were unexpected. Most especially, we have noted 
changes in the geometry of the circumstellar environ- 
ment as reflected in the disappearance of one of the two 
columns, and affecting ~ 25% of the hard emission. The 
strong attenuation of the 2001 soft X-ray spectrum was 
clearly evident. In addition, the strengths of the K fluo- 
rescence emission features decreased in line with the de- 
crease of the column absorption, suggesting that the part 
of the emission is caused by scattering of hard photons 
through the attenuating column that was present at the 
earlier time. 
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The reduced attenuation of the soft X-rays partially 
accounts for the improved ability to study the nature of 
the warm and cool plasma emissions. As a result, it is fi- 
nally clear that these components are almost monothermal 
and therefore are not part of an integral structure with a 
smoothly varying thermal emission measure, such as an 
accretion column or a cooling flow plasma. Moreover, the 
broadening of the lines has increased during the 2001- 
2004 interval, cither because of an increase in a quasi- 
turbulence or the splitting of a former single region into 
two with different projected radial velocities. 

In contrast to the warm component, a study of the vari- 
ations in the Fe K-shell lines with temperature shows that 
the hot plasma need not consist of a uniform fcTg value. 
Rather, we believe that small variations in the hardness 
that are occasionally observed in our data and those dis- 
cussed by SRC98, RSOO, and RSH02 can be understood 
as variations of the instantaneous average temperature re- 
sulting from the rapid evolution of small number of flares 
and nearby basal emission regions. 

The unique properties of each of the thermal com- 
ponents, the changing absorption column geometry, and 
increases in line broadening all provide new hints to 
the mechanism responsible for the X-ray production. 
Although on one hand, the disappearance of the strong 
absorption column and the lack of correlation between the 
fluorescence features no longer support the argument that 
the gas in the Be star's circumstellar disk strongly inter- 
act with the hard X-rays, the discreteness of the plasma 
components argues that they are likely to occupy distinct 
volumes. 

Future generations of X-ray telescopes such as the 
International X-ray Observatory (IXO) will be important 
in refining our understanding of the spatial distribution 
of emission volumes of 7 Cas and its analogs. Their obser- 
vations promise to address such question as whether the 
observed line broadening can be tied to the rotational ve- 
locity of the Be star and to resolving distinct sources, such 
as corotating active regions. 
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